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Composite dark matter candidates, which can arise from new strongly-coupled sectors, are well- 
motivated and phenomenologically interesting, particularly in the context of asymmetric generation of 
the relic density. In this work, we employ lattice calculations to study the electromagnetic form factors of 
electroweak- neutral dark-matter baryons for a three-color, QCD-like theory with Nf — 2 and 6 degenerate 
fermions in the fundamental representation. We calculate the (connected) charge radius and anomalous 
magnetic moment, both of which can play a significant role for direct detection of composite dark mat- 
ter. We find minimal Nf dependence in these quantities. We generate mass-dependent cross-sections for 
dark matter-nucleon interactions and use them in conjunction with experimental results from XENON100, 
excluding dark matter candidates of this type with masses below 10 TeV. 



PACS numbers: ll.10.Hi, ll.15.Ha, 95.35.+d 

Introduction Experimental bounds on the interaction 
of the dark matter with Standard-Model (SM) particles 
have strengthened by many orders of magnitude in recent 
years. In particular, dark-matter particles cannot have SM- 
strength couplings to electroweak gauge bosons, based on 
direct-detection constraints HIE). At the same time, there 
is a strong motivation for the dark matter to couple to the 
SM in some way for the purpose of relic density generation, 
either as a thermal relic via the so-called "WIMP miracle" 
(see J51 for a recent review) or through an asymmetric sce- 
nario which may be related to the creation of baryon asym- 
metry iHHU- Construction of dark matter models thus re- 
quires a careful balance between the presence and absence 
of dark-sector interactions with the SM. 

Composite dark matter models provide a simple mech- 
anism for attaining this balance, one which can lead to 
interesting and unique phenomenology. By hypothesiz- 
ing a new, confining gauge force in the dark sector, an 
electroweak-neutral composite dark matter candidate can 
be constructed as a bound state of electroweak-charged 
constituents. In this way, electroweak interactions can be 
active in the early Universe for the generation of relic den- 
sity, but only neutral bound states survive to the present 
day. Electroweak coupling to the constituents is still possi- 
ble, leading to form-factor suppressed interactions with the 
neutral composites. They can be roughly estimated from 
QCD analogs, but in general can be determined quantita- 
tively only by lattice calculations. 

In this paper, we consider an underlying SU(3) gauge 
theory with fermions in the fundamental representation, but 



focus on fermions not associated with electroweak break- 
ing. We use SU (3) because much is known about it from 
lattice QCD and because we have already generated lattice 
vacuum states of SU(3) with 2 and 6 fundamental flavors 
on large lattices ll9l [T0ll . We take the fermions to be mass- 
degenerate SU (2) l singlets such that Q = Y. We con- 
sider a two-fermion theory (Nf = 2) with Q u = 2/3 and 
Q d = —1/3, as well as a six-fermion theory (Nf = 6) 
with three such pairs of fermions. In either case, the light- 
est baryon is expected to be electrically neutral, and will 
therefore also have vanishing weak charge. The dominant 
contribution to its interaction with ordinary nuclei will be 
due to single photon exchange, which can be parameter- 
ized primarily in terms of its magnetic moment and charge 
radius. In these initial lattice calculations we consider only 
quark-line connected contributions to the charge radius and 
magnetic moment. We compute the electromagnetic form 
factors of this particle to extract these quantities, describe 
their dependence on Nf, and discuss consequences for di- 
rect detection. 

One could also modify or enlarge the fermion content 
of the SU(3) gauge theory to include S'C r (2) £ -doublet 
fermions. This would be a necessary modification in or- 
der to consider composite dark matter arising in a theory of 
dynamical electroweak symmetry breaking PM6l fTTT - iT9l . 
Careful model building is then required to ensure that the 
lightest baryon is net electroweak neutral. We do not dis- 
cuss this possibility here. 

Model setup For the theory with 5'C/(2)i-singlet 
fermions carrying charges Q u = 2/3 and Q d = —1/3, 
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with Nf = 2 or 6, the analogue of the neutron (N ~ udd) 
will be the dark matter candidate, with mass M b and car- 
rying no net electroweak charge. It is stabilized by conser- 
vation of dark baryon number. The other charged baryons 
are expected to be heavier due to electromagnetic mass cor- 
rections of order AM ~ aM s /47r. We include a fermion 
mass m f, essential for lattice calculation purposes, and ex- 
amine dependence on m/ for a range rrif M B - 

Our dark sector also contains Nj — 1 pseudo-Nambu- 
Goldstone -boson (PNGB) states. We assume that these 
states are unstable, decaying to Standard-Model particles 
with a sufficient rate that their presence does not influence 
the cosmological history of the Universe. 

As our focus is on direct-detection signatures, we do not 
consider the dark matter generation in detail here. The con- 
finement scale A, or equivalently the dark matter mass M B , 
is a free parameter in our construction. 

Electromagnetic Form Factors Since the neutral 
baryon in the SU (2) -singlet theory is the dark matter can- 
didate of interest |[32l . the baryon mass M b (degenerate in 
the absence of other interactions) is the dark matter mass. 
This mass and all other dimensionful quantities are ex- 
pressed in lattice units here. 

The quantities of central interest here are the Dirac and 
the Pauli electromagnetic form factors of a neutral dark- 
matter baryon \N(p)). For the Nf = 2 case, they can be 
expressed in terms of matrix elements of the vector currents 
of individual quarks as follows: 
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where ip = u, d are quark fields, U, U are on-shell baryon 
spinors, q = p' — p, and Q 2 = —q 2 > is the mo- 
mentum transfer. In the forward limit Q 2 = 0, the Dirac 
form factors are equal to the numbers of the valence quarks: 
F?(0) = landF^O) = 2. 

From these one constructs the isovector and isoscalar 
form factors ll33l : 
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Both of these quantities can be extracted from lattice cal- 
culations, but the isoscalar contribution contains expensive 
disconnected lattice quark contractions, which cancel in the 
isovector case, and as a result, isovector form factors are far 
more tractable. While we ultimately will calculate the dis- 
connected pieces of the isoscalar form factor as well, this 
work will focus on only the connected contributions. 

For the Nf = 6 case, with three pairs of u(Q = 2/3) 
and d(Q = —1/3) fermions, we take the \N(p)) state 
to be composed of fermions from only one pair. Since we 
omit disconnected lattice quark contractions in our calcula- 
tion, it is only the currents ip^^ip composed of the fermion 



fields from the same pair that contribute to the computed 
electromagnetic form factors. Therefore, in our calculation 
the other two pairs play a role in only the strong dynamics 
of the SU(3) gauge theory. 

The full electromagnetic form factors of the neutral dark 
barvon ll34l are given by 

-Fl,2;neut(Q 2 ) = Q u Fi 2 {Q 2 ) + QdFi 2 {Q 2 ) 
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since F?(0) = 3 and Ff(Q) = 1, the total charge 
-Fi;neut(0) = 0. For soft single -photon exchange scatter- 
ing, only the forward (Q 2 — > 0) behavior of the electro- 
magnetic form factors is relevant. Since the electric charge 
-F\;neut(0) is zero, only the magnetic moment ^ neu t = K neu t 
and the Dirac radius (r 2 . neut ) contribute to the scattering 
amplitude to the lowest order in Q 2 : 

1 
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t + 0(Q 2 ), 
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The Dirac charge radius (r^. neut ) determines the slope of 



the form factor in the Q 2 
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> limit: 
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The definition of the radius ([5]) is motivated by the alge 
braic identity 

,3 2 / \ _ a dF!(Q 2 
d rr p(r) = — o- 



dQ 2 
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where p(r) is the "charge density" 



d 3 re iq ~p(r) = F^Q 2 ) , Q 2 
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which has physical meaning if and only if the spatial extent 
of this distribution is much larger than the Compton wave 
length of the composite particle, (r 2 ) S> Mg 2 . Since the 
total charge, J d 3 r p(r) = -Fi(O), is zero, the charge den- 
sity must have alternating sign (or be exactly zero), and the 
integral in Eq. (|6| can be either positive or negative. 

For the following, we also need to define the mean 
squared charge radius {r 2 E ), or the "radius" of the charge 
form factor G E {Q 2 ), 



G E (Q 2 ) = F 1 (Q 2 
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Similar to Eq. ^5), the charge radius of the neutral baryon 
is equal to 

dG £ ; neut(Q 2 ) 

dQ 2 Q2=0 
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differing from the Dirac radius by only the relativistic cor- 
rection ~ Mq (the Foldy term). This correction is impor- 
tant if the size of the particle is comparable to its Compton 
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wave length, which is the case for the neutron and the pro- 
ton in QCD. 

The (anomalous) magnetic moment of the neutral baryon 
is related to the isovector and isoscalar moments as 

1 1 (.0) 
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The isovector and isoscalar Dirac form factors are not zero 
in the forward limit. Their radii are defined to be indepen- 
dent of their overall normalization, 
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The radii of the neutral baryon are related to the isovector 
and isoscalar radii as follows: 

(rLeut) = \{rlY ~ \{r\Y , 

(12) 
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Simulation Details Lattice calculations are performed 
using 32 3 x 64 domain-wall lattices with the Iwasaki im- 
proved gauge action and a fifth-dimensional length L s = 
16 and a domain-wall height of m Q = 1.8. By us- 
ing domain-wall fermions, the calculation preserves ex- 
act flavor symmetry, and chiral-breaking lattice spacing ar- 
tifacts are suppressed. The calculation is performed for 
N f = 2 at p = 2.70 and N f = 6 at = 2.10. 
The beta values are tuned to match the confinement scale 
of both theories relative to the lattice spacing, includ- 
ing M B as we shall see below. For both Nf = 2 and 
Nf = 6, five separate mass points are analyzed with mt = 
0.010,0.015,0.020,0.025,0.030. The pion masses (in 
units of the nucleon mass) are 0.41 < m^jM B < 0.52 
and 0.44 < m„/M B < 0.52 for N f = 2 and N f = 6, 
respectively. Further details and other results from these 
ensembles are given in ll9l [T0ll20l . 

Calculation and Fitting The parameters of interest are 
extracted from two sets of correlation functions: two-point 
correlation functions given by 



VA r(r,p) = ^e-' ip - x (iV(x,r)iV(0)), 
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and three-point correlation functions 
C NO n(t, T,p,p' 
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where 0(y, r) is the quark vector current density operator. 

The long-distance limit of the Euclidean time behavior 
of these correlation functions is given by 
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where r po i is the polarization matrix of the initial and final 
baryon spin states corresponding to Eq. ( 13|14 i, T M is the 
fermion vertex function (cf. Eq.Q), 

r» = F.mr + F 2 m^ , (i7) 

and A is the difference in energy between the ground 
and the first excited state of the baryon. More details on 
the form factor calculation on the lattce can be found in 
Ref. ED- 

It is useful to form an appropriate ratio of these correla- 
tion functions 
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where the long Euclidean time behavior yields 
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In general, the excited state corrections can lead to signifi- 
cant systematic errors on three-point functions [(22). 

The form factors Fi^(Q 2 ) are calculated at discrete val- 
ues of the momentum transfer Q 2 ~ (p' — p) 2 determined 
by the lattice volume. We interpolate the Dirac and isovec- 
tor Pauli form factors using a dipole formula fit 



^1,2 (Q 2 
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motivated by nucleon form factor phenomenology. The 
isoscalar Pauli form factor turns out to be very close to 
zero, and the dipole form that has definite sign does not 
necessarily yield a stable fit to the data; therefore, we use 
the linear fit F^Q 2 ) ~ F J (0)+ J F 1 /(0)Q 2 . We use these fits 
to interpolate (extrapolate in the case of Pauli form factors) 
near the forward limit Q 2 = in order to determine k and 

<r?>. 
Lattice Results 

Baryon Mass The dark-matter baryon mass is plotted 
as a function of the fermion mass nif in Fig. [T] A lin- 
ear dependence of the baryon mass on mf can be seen for 
both theories, as expected in the calculation regime where 
the fermion masses are small. In the absence of additional 
interactions, a finite value of my is required to give mass 
to the PNGB's of the theory, but we nevertheless perform 
a linear fit in order to extract the chiral-limit baryon mass 
M Bo . This scale, which can be taken as a proxy for the 
confinement scale of the theory, serves as a common refer- 
ence scale for the calculation results with nif > 0. 
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FIG. 1 : Dark-matter baryon mass (in lattice units) with Nf = 2 
(red) and Nf — 6 (blue), as a function of the fermion mass m,f 
(also in lattice units). The two data sets are extrapolated to obtain 
the chiral-limit baryon mass M Bo , which is used to set a physical 
scale independent of arrif. With the chosen lattice couplings, 
M B(j is the same within statistical precision in the Nf — 2 and 
Nf — 6 theories. 
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FIG. 2: The neutral baryon anomalous magnetic moment for 
Nf — 2 (red) and Nf — 6 (blue) theories versus dark-baryon 
mass. This quantity shows no systematic separation between two 
and six flavor theories. 



Anomalous magnetic moment The anomalous mag- 
netic moment is the most important for direct detection ex- 
periments. It enters at the dimension-5 level in the baryon 
effective field theory and arises as the zero-momentum 
value of the Pauli form factor, ^(0). The iso vector Pauli 
form factor, giving k v , is under most control since all ex- 
pensive disconnected contributions cancel due to isospin 
symmetry. The isoscalar channel, which is also necessary 
to determine K n eut, has both connected and disconnected 
contributions to the three-point correlation function. In this 
initial work, we omit the disconnected contributions and 
assume the connected pieces dominate the isoscalar contri- 
bution as observed in lattice QCD. 

We plot the anomalous magnetic moment n Deut , com- 
puted as described above, versus M B /M Bo in Fig. [2] It 



shows little dependence on the mass and little dependence 
on the number of fermions. The Nf = 2 results K neut ~ 
— (1.71 . . . 2.09) are consistent with the measured neutron 
value k = —1.91 ||23ll . Calculations of nucleon structure 
with Nf = 2 Wilson fermions were previously reported in 
Ref.[24|, which found values K neut ~ —(1.30 . . . 1.45). 
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FIG. 3: The neutral baryon mean squared charge radius (in lattice 
units) for Nf = 2 (red) and 6 (blue), versus dark-baryon mass. 
Again, no significant systematic difference between the two the- 
ories is seen over the range of masses considered. 

Charge radius While the charge radius is expected to 
lead to a smaller effect on the spin-independent cross sec- 
tion as compared to the magnetic moment, it could have a 
significant effect if its value depends significantly on Nf. 
It is therefore informative to explore the relative size of the 
charge radius contribution to the spin-independent cross 
section. As with the magnetic moment, only the isovec- 
tor charge radius is absent of disconnected lattice quark 
contractions, but we omit them for the isoscalar channel 
as well. 

The results for the mean square Dirac charge radius 
i r E -neut) °f an electroweak-neutral dark-matter baryon are 
presented in Fig. [3] Note that the results are negative (see 
discussion after Eq. Q). As in the case of the anomalous 
moment, our results show little dependence on Nf and lit- 
tle dependence on the dark-baryon mass as it varies due 
to changes in the underlying fermion mass. If the fermion 
mass is reduced further, bringing M B /M B(i closer to unity, 
the magnitude (r|;. neut ) is expected to grow. This is be- 
cause the PNGB mass drops, and the charge radius is quite 
sensitive to the size of the PNGB cloud. 

For Nf = 2, this point can be made more pre- 
cisely by comparison to QCD. There, the mean squared 
charge radius of the neutron is also negative, (r% n ) = 



-0.1161 22 fm z 



Our N 



f 



2 calculation cor- 



responds to QCD with M B ~ 1 GeV and lattice spac- 
ing a ~ 0.055 fm, but with relatively heavy underlying 
quarks, and thus relatively heavy pions: the pion mass in 
units of M B ranges between the lightest m^/m B = 0.41 
to the heaviest m„/m B = 0.52. In QCD units, our re- 
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suit is (r| neut ) « -(0.009 . . . 0.025) fm 2 , substantially 
less than the observed result. Previous calculations of nu- 
cleon structure with Nf = 2 Wilson fermions 11241 yielded 
similar values (r| neut ) = -(0.011 . . . 0.023) fm 2 . These 
results, too, employed relatively heavy underlying quarks. 
In our case, further studies with smaller fermion mass can 
shed light on the range of direct detection allowed values 
for the mean square charge radius. 
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FIG. 4: Calculated XENON100 event rates based on energy cuts 
and acceptance rates from Ref. 1251 (solid lines). For comparison, 
we also show scattering rates using only the charge radius term, 
which is suppressed by two additional powers of M x (dashed 
lines). The experimental upper bound on event rates, based on 
accumulated 2323.7 kg-days of exposure 1 25 1 are shown with the 
dotted lines. 



Direct detection exclusion plots We next compare our 
calculations of dark-matter parameters with the current 
experimental bounds on the dark-matter-nucleus cross- 
sections in direct detection experiments. Currently, the 
most stringent bound is provided by the XENON100 ex- 
periment 11251 . in which hypothetical dark-matter particles 
are detected through their collisions with xenon nuclei with 
Z = 54 and A = 124 . . . 136, and which has accumulated 
2323.7 kg-days of effective exposure. Two of the isotopes, 
129 Xe and 131 Xe, have non-zero spin and are sensitive to 
the spin-dependent Ml interaction. Their combined abun- 
dance constitutes approximately 1/2 in natural xenon 11261 . 

In this section, we adopt a more conventional notation 
M x for the mass of the dark-matter particle, and also de- 
note its radius and magnetic moment with a subscript "x". 
Figs. [2] and [3] show that the anomalous moment and 
mean square charge radius vary little with the amount of 
the dark-matter mass coming from the underlying fermion 
mass (and also vary little as Nf is increased from 2 to 6). 

The differential cross-section of a dark-matter fermion 
and a nucleus, to leading order in the non-relativistic dark- 



matter velocity v <C 1 is 
da 



|-Msi| 2 + |-Msd| 



dE R 16vr(M x + M T ) 2 E™ 



(21) 



where M T is the mass of the target nucleus, and E^ = 
(m x +Mt) 2 * s trie max i ma l recoil energy for given colli- 

X 2 

sion velocity v. The quantities and |A^si,sd| are spin- 
independent amplitudes squared, averaged over initial 
and summed over final states: 
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Here, Z and A are the charge and atomic numbers of 
a specific xenon isotope, (ht/ fJ, n ) is the nucleus magnetic 
moment expressed in Bohr magnetons \i n = , F C ^ S (Q) 
are its nuclear charge and spin form factors, respectively, at 
the momentum transfer Q ~ \/Q 2 = \/2M T £ , fl , and #cm 
is the scattering angle in the center-of-mass frame ll35l . For 

non-relativistic velocities, cot 2 = 



1 



For the nuclear form factors F CjS {Q 2 ), we use the fol- 
lowing phenomenological expressions lfT7ll27l : 



l^ c (Q)| 2 = 9 



sm(QR c ) — (QR C ) cos(QR c 



\Fs(Q)\' 



0.047, 

sin(Qfl s ) 



{QRcY 



2.55 < QR S < 4.5, 



<Qsy 



(24) 



, otherwise , 
1.14^ 1 / 3 fm, R s = 1.00 A 1 ' 5 fm, and S 



(25) 



where R c 
0.9 fm. 

Following Refs. lfT71l25l , we compute the scattering rate 
for a range of dark-matter particle masses with the recoil 
energies E R = 6.6 . . . 43.3 keV: 

R = — / dE R Acc(E R ){< 



da 



dE R /f J 
(26) 

where (•) f denotes averaging over the DM velocity distri- 
bution (27 1, v' = \v — ^Earthl is the dark-matter velocity 



with respect to the detector (the Earth), and Acc(E R ) is 
the recoil energy-dependent acceptance rate of the detec- 
tor ll25l . We assume the thermal distribution of velocities 
in the galactic dark-matter halo 11271 . 

(27) 

with v = fEarth = 220 km/s, v esc = 544 km/s, and the 
dark-matter mass density /3dm — 0.3 GeV/cm 3 . Finally, 
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we average the expected scattering rate over the natural 
xenon isotopic abundances. 

We show computed scattering rates in Fig. H] with solid 
lines. The accumulated XENON100 statistics 1251 exclude 
composite dark matter particles with M x < 10 TeV. With 
fixed values of the dimensionless quantities M 2 (r 2 ) and k 
computed on a lattice, the differential cross-section scales 
as 



da 



A 



M 2 



(28) 



The charge radius contribution is suppressed as M~ 2 rel- 
ative to the magnetic moment contribution and becomes 
negligible with growing M x . In the scattering rate shown 
in Fig. [4] both contributions are additionally suppressed by 
the DM particle number density Pdm/M x as M x ~~ ^ 00 
(see Eq.([26]l). The large-M x scaling of the charge radius 
term is sEown in Fig. [4] with the dashed lines; it is ev- 
ident that the total scattering rate (solid lines) is domi- 
nated by the magnetic moment term for dark matter masses 
M x > 25 GeV. Even if one were to make the charge ra- 
dius as much as an order of magnitude larger by reducing 
the PNGB mass (see the discussion following Fig. |3j, its 
contribution would still be negligible at M B = 10 TeV, 
the lower limit of the allowed region. 

Discussion We have studied the electromagnetic form 
factors of electroweak-neutral dark-matter baryons in an 
SU(3) gauge theory with N f = 2 and 6 S*7(2) L -singlet 
fermions, with charge assignments +2/3 and —1/3 (one 
pair or three pairs). These baryons have the desired proper- 
ties of dark matter since they are stable, electroweak neu- 
tral, and can explain the relic density through the same 
early universe sphaleron process that describes baryogene- 
sis. 

Of particular interest to direct detection experiments are 
the anomalous magnetic moment and mean square charge 
radius of the dark-matter baryon. These parameters deter- 
mine the observed cross-section with nuclei (in this work, 
we primarily focus on Xenon) due to (dominant) single- 
photon exchange. The contribution from the dark-matter 
anomalous moment dominates the charge-radius contribu- 
tion for M v > 25GeV. However, in our calculation the 
charge radius {r 2 E ) turns out to be particularly small, much 
smaller than it is in QCD. Exploring smaller quark mass re- 
gions may change this balance and make the charge radius 
more relevant for the direct detection of the dark matter. 

Examining the dark matter exclusion plots in light of the 
latest Xenon 100 results fl25l , we conclude that in these the- 
ories, dark-matter masses less than 10 TeV are excluded. 
We have so far seen little dependence on Nr. It will 
be interesting to see whether this begins to change, even 
continuing to neglect disconnected quark contractions, as 
Nf is increased toward the edge of the conformal window 
(Nf ~ 10 — 12 for an SU(3) gauge theory with fermions 
in the fundamental representation). When the disconnected 



contractions are included, additional Nf dependence will 
arise simply from the counting of these loop contributions. 

As we have shown in this work, with non-zero magnetic 
moment the experimental constraints on the dark matter 
mass are quite stringent. This naturally motivates the con- 
sideration of even N c theories, in which the baryons are 
bosonic and thus have no magnetic moment. Interactions 
can be further suppressed if the charge assignments are 
symmetric in such a way that the charge radius vanishes 
(see e.g. |[T8l ). making the electromagnetic polarizabilities 
the dominant interactions. 
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